clear all
load W_time1_age0_eh1_comp11.out
n = 50;
w_f_g=W_time1_age0_eh1_comp11; % Read wages, g and f distributions
wf_1=w_f_g(:,1);        % Wage in formal sector
wi_1=w_f_g(:,2);        % Wage in informal sector
wf_2=w_f_g(:,7);      % Wage in formal sector
wi_2=w_f_g(:,8);      % Wage in informal sector

load Rw_1_f_gamma_0_GE Rw_1_f_gamma_0_GE
load Rw_2_f_gamma_0_GE Rw_2_f_gamma_0_GE
load Rw_1_i_gamma_0_GE Rw_1_i_gamma_0_GE
load Rw_2_i_gamma_0_GE Rw_2_i_gamma_0_GE

% Individual model (women)
load Rw_2_f_ind
load Rw_2_i_ind

x_head = wi_1; x_spouse = wi_2;

plot(x_head,Rw_2_i_gamma_0_GE,'-b',x_head,Rw_2_f_gamma_0_GE,'-g',x_head,Rw_2_i_ind*ones(50,1),'--b',x_head,Rw_2_f_ind*ones(50,1),'--g','LineWidth',2.5)
 
xlabel('Wage of men','FontSize',15)
ylabel('Reservation wage of women','FontSize',15)
box off
 
legend('$\hat{w}_{in->ii}$','$\hat{w}_{in->if}$','$\hat{w}_{n->i}$','$\hat{w}_{n->f}$','Interpreter','latex','FontSize',20)


print -dpdf Rw_ind_hh_search_women.pdf

